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Abstract. Outflows from quasars inevitably pollute the intergalactic medium (IGM) 
with magnetic fields. The short-lived activity of a quasar leaves behind an expanding 
magnetized bubble in the IGM. We model the expansion of the remnant quasar bub- 
bles and calculate their distribution as a function magnetic field strength at different 
redshifts. We find that by a redshift 2 ~ 3, about 5-80% of the IGM volume is filled by 
magnetic fields with an energy density > 10% of the mean thermal energy density of a 
photo- ionized IGM (at ~ 10* K). As massive galaxies and X-ray clusters condense out 
of the magnetized IGM, the adiabatic compression of the magnetic field could result in 
the fields observed in these systems without a need for further dynamo amplification. 

1 Introduction 

Clusters of galaxies contain substantial magnetic fields with strengths B ^ 0.1 — 
10 /iG and coherence lengths ^ '--^ 10 kpc 0]. The origin of such fields could have 
important implications for structure formation. Assuming flux conservation, a 
cluster field Bc\ ~ 10~^ G would imply Bigm ^ 10^^ G in the diffuse intergalac- 
tic medium (IGM), which would constitute ~ 5% of the thermal energy density 
of a photoionized IGM. The observational constraints on an intergalactic mag- 
netic field (IGMF) are weak, requiring only that Bjgm < 10"^(^/ Mpc)-^/^ G 
in the currently popular ylCDM model [Q. 

Unfortunately, there is no convincing model for the formation of cluster fields 
or the IGMF. Because the rotation times of clusters exceed the Hubble time, 
dynamos are expected to be ineffective. Primordial generation scenarios cannot 
explain the large coherence lengths |Q. Models that generate magnetic fields 
during structure formation or in starbursts Q suffer from similar problems. 

We have examined the possibility that the IGMF was originally produced 
near supermassive black holes and expelled into the IGM through mechanical 
outflows 1^ from radio-loud quasars (RLQs) and broad absorption line quasars 
(BALQs). Supermassive black holes are one of the few classes of astrophysical 
objects with energy reservoirs large enough to account for the large-scale fields 
in clusters, and the relatively small number of powerful sources can accomodate 
the large observed coherence lengths, as emphasized by . For more details on 
the model, see 0. 
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2 Filling Factor of Magnetized Regions &; Mean Fields 

While a quasar is active, its outflow is in the form of twin coUimated jets (for 
a RLQ) or equatorial winds (for a BALQ). After the quasar becomes dormant, 
the outflow remnant is overpressured with respect to the IGM and continues 
to expand in comoving coordinates until its outward velocity matches the Hub- 
ble flow velocity. We assume that the outflow remnant isotropizes and expands 
adiabatically as a spherical shell during this late phase. 

Simple energy conservation provides a surprisingly accurate estimate of the 
final comoving bubble size i?max- Balancing the energy input of the quasar and 
the final kinetic energy of the shell, we find that i?max oc [LqTqEKi^ + £b)Y^^, 
where Lq is the luminosity of the quasar, Tq is its lifetime, Ek is the ratio of the 
mechanical and radiative luminosities, and Eb is the ratio of the magnetic and 
mechanical energy outputs. Note the weak dependence on the quasar parameters, 
making our results robust to uncertainties in their measurement. In particular, 
we find that magnetic fields do not play an important role in the expansion for 
realistic values of Eb- We therefore ignore the geometric and magnetohydrody- 
namic effects of the magnetic field and assume simple flux conservation. 

We next calculate the number of quasar sources. For z < 4, we use the ob- 
served optical luminosity function of quasars |^ together with an assumed inci- 
dence rate of outflows /. For z > 4, we assume that the incidence rate of quasars 
is proportional to the Press-Schechter mass function, with the proportionality 
constant set by matching to the observed luminosity function at z ^ 4 |9| . 

We examine two quasar models: the RLQ model, with eb ~ 0.1 and / = 
0.1 ||l§, and the BAL model, with eb = 0.01 and / = 1 Q. Note that there 
are no existing observations of magnetic fields in BALQs, though the fields may 
play an important role in the wind mechanism p^ . 

We then calculate the total filling factor of the quasar bubbles as a function of 
redshift, F{z), and the global volume-averaged magnetic energy density, ub{z), 
for our two models by numerically integrating the equation of motion of each 
remnant and summing over all quasar sources. Our results are shown in Fig. 1, 
along with the distribution of magnetic field strengths for a series of redshifts. 



3 Results 

We predict a cellular IGMF filling a substantial fraction of space by z ~ 3 (Fig. 
1), with each cell a fossil bubble produced by a single magnetized quasar outflow. 
Cells generated by RLQ quasars are more highly magnetized than those from 
BAL quasars but fill a smaller volume. We predict that i?iGM ~ 10/e_B nG when 
averaged over all of space; for each of our models, simple adiabatic compression 
of such an IGMF can account for the observed cluster magnetic fields. 

Direct detection of the cells via Faraday rotation measurements will be diffi- 
cult, but electron acceleration by shocks in the cells will cause synchrotron emis- 
sion. The same accelerated electrons produce 7-rays through inverse-Compton 
scattering of the cosmic microwave background. Correlated maps of the 7-ray 
and radio skies may allow us to calibrate the magnetic field in the shocks JTsf. 
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Finally, we note that the non-thermal pressure of the magnetic fields may 
help to resolve the discrepancy between simulated and observed line widths in 
the Lya forest Q. 
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Fig. 1. Top Left: Volume filling fraction of magnetized bubbles F{z). Bottom Left: 
Ratio of magnetic energy density, ub, to the fiducial thermal energy density Ufid = 
3n{z)kTiGM , where Tigm = 10'* K. In each of these panels, results are shown for both 
the RLQ and BALQ models. Right: Probability distributions of bub ble mag netic field, 
P{B,z), for the RLQ model, at various redshifts. Here B = B/{y^eB /O.l); in these 
units the curves are independent of es. The dotted curve shows the contribution to 
the z = distribution function from quasars forming at 2:0 < 4 and the dashed curve 
shows the contribution from quasars at zo > 4. All curves are normalized to have unit 
area. 



